Hidden Reentrant and Larkin-Ovchinnikov-Fulde-Ferrell Superconducting Phases in a 

Magnetic Field in (TMTSF)2C104 
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We solve a long-standing problem about a theoretical description of the upper critical magnetic 
field, parallel to conducting layers and perpendicular to conducting chains, in (TMTSF)2C104 su- 
perconductor. In particular, we explain why the experimental upper critical field, H^2 — 6 T, is 
higher than both the quasi-classical upper critical field and Clogston paramagnetic limit. We show 
that this property is due to the coexistence of the hidden Reentrant and Larkin-Ovchinnikov-Fulde- 
Ferrell phases in a magnetic field in a form of three plane waves with non-zero momenta of the 
Cooper pairs. Our results are in good qualitative and quantitative agreement with the recent ex- 
perimental measurements of H^2 S'lid support a singlet d-wave like scenario of superconductivity in 
(TMTSF)2C104. 
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Physical properties of quasi-one-dimensional (QID) or- 
ganic conductors (TMTSF)2X (X^PFg, CIO4, Re04, 
etc.) have been intensively studied [1,2] since a discovery 
of superconductivity in (TMTSF)2PF6 [3]. Early exper- 
iments [4,5] clear showed that superconducting phases 
in these compounds were unconventional and that the 
corresponding order parameters changed their signs on 
QID Fermi surfaces (FS). In particular, it was shown 
that the Hebel-Slichter peak was absent in the NMR 
experiment [4] and superconductivity was destroyed by 
non- magnetic impurities [5]. These results have been re- 
cently confirmed in a number of publications (see, for 
example, Refs. [6,7]). The first Knight shift measure- 
ments [7,8], performed in (TMTSF)2PF6 conductor in a 
magnetic field H = 1.43 T, showed that the Knight shift 
was unchanged in superconducting phase and were in- 
terpreted as evidence for triplet superconductivity. On 
the other hand, more recent Knight shift data [9], per- 
formed in (TMTSF)2C104 conductor, clear demonstrate 
the Knight shift change through the superconducting 
transition in a magnetic field H = 0.957 T. They are 
interpreted [9] in terms of singlet pairing in superconduc- 
tor (TMTSF)2C104 at least at relatively weak magnetic 
fields. 

Another source of information about a spin part of 
the superconducting order parameter was provided by 
fact that the experimental upper critical magnetic field 
along conducting chains, iJ°2 [10], was clear paramag- 
netically limited [11]. This has been recently confirmed 
in Refs. [12-14]. In addition, new superconducting phase 
has been discovered in (TMTSF)2C104 [12,13] for a mag- 
netic field, parallel to conducting chains. The suggested 
hypothesis [12,13] that it can be the Larkin-Ovchinnikov- 
Fulde-Ferrell (LOFF) phase [15,16] has been recently the- 
oretically supported [17]. Note that the above men- 
tioned experimental and theoretical works are in favor 
of a singlet d-wave like scenario of superconductivity in 
(TMTST)2C104 [14,17-19] . 



In this situation, where support for a singlet 
d-wave like scenario of superconducting pairing in 
(TMTST)2C104 conductor is growing, it is important 
theoretically reinvestigate the upper critical field, paral- 
lel to conducting layers and perpendicular to conducting 
chains, -ff^2- Fo'^ many years, large experimental values 
of H^2 [10,12,13,20-23], which exceeds both the quasi- 
classical upper critical field, i?c2(0) [24], and Clogston 
paramagnetic limit. Hp [25], have been considered as one 
of the main arguments in favor of triplet superconductiv- 
ity. Although the exceeding of the values of i?c2(0) and 
Hp was predicted for H^2 in both singlet and triplet cases 
[26-30], it was also shown that, for realistic band param- 
eters of (TMTST)2X conductors, it can happen only in 
a triplet case [11,20,26-30]. 

The goal of our Letter is to demonstrate that super- 
conductivity in (TMTSF)2C104 can exceed both criti- 
cal magnetic fields, H^2'S^) ~ 3.5 T and Hp ~ 2.7 T, 
and reach its experimental value, H^2 — 6 T [12,13,23], 
even in case of a singlet d-wave like superconducting 
pairing. The first our point is that the Pauli param- 
agnetic effects in all previous theories [26-30,11,20] were 
treated not completely correctly. The second our point 
is that the 3D — 2D dimensional crossover [26] hap- 
pens at magnetic fields H^ ~ 5 — 6 T, which are much 
lower than the previously assumed. The latter state- 
ment is shown to result from theoretical analysis of both 
the Ginzburg-Landau (GL) slopes, dH^2/'^'^\T=T^ and 



dH^2/ dT\T=Tc^ measured in Refs. [12,13,23], and the so- 
called Lee-Naughton-Lebed oscillations [31,32]. In the 
Letter, we derive a novel gap equation, which treats ac- 
curately both the Pauli paramagnetic and orbital de- 
structive effects against superconductivity. By analyzing 
this equation, we show that it predicts the upper criti- 
cal field, H^2 — 6 T, for real values of band parameters 
in (TMTSF)2C104. Superconducting phase, which ex- 
ists at such high magnetic fields, is shown to be very 
peculiar. It is characterized by a inhomogeneous order 
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parameter in a form of the three LOFF-like waves, which 
appear both due to the 3D — ^ 2D dimensional crossover 
and PauU paramagnetic effects. It is important that this 
phase is characterized by the Cooper pairs, localized on 
conducting layers, with probability of the Cooper pair 
jumping from one layer to another being small. There- 
fore, it is not destroyed by the orbital effects in a parallel 
magnetic field. In the absence of the Pauli paramagnetic 
effects, such phase would correspond to the Reentrant 
superconductivity with dTc/ dH > 0, therefore, we call it 
the hidden Reentrant superconducting phase. 

Below, we consider a tight-binding orthorhom- 
bic model of anisotropic QID electron spectrum in 
(TMTSF)2C104 conductor, 

e(p) = -2taCos{pxa/2) - 2tbCos{pyb') - 2icCos(p^c*), 

(1) 

which can be simplified near two slightly corrugated 
sheets of QID FS as 

^e'''(p) = ±'Vx{Pv)\Px T Pf{Pv)] - ^tc cos{pzC*). (2) 

[Here ta ^ U ^ tc correspond to electron hoping inte- 
grals along a , b', and c* axes, respectively; -|-(-) stands 
for right (left) sheet of FS.] 

In a magnetic field, parallel to conducting planes and 
perpendicular to conducting chains of QID conductor. 



H=(0,i7,0), A = (0,0,-7Jx), 



(3) 



we use the so-called Peierls substitution method, Px =F 
ppiPy) — >■ —id/dx, Pz ^ Pz — &Az/c. As a result, effec- 
tive Scrodinger equation for electron wave functions in a 
mixed representation, 'tjj^{x,Py,Pz,cr), can be written as 



-2tcC0s[ PzC* H X —iiBoH 

Vf 



X ip^ {x,py,pz,(T) = Se ipf {px ,Py,Pz, i^) 
with electron wave in a real space functions being 
"^tix, y, z, a) = ex.p[ipF{Py)x] exp{ipyy) ex.p{ipzz) 



i)f{x, 



(5) 



where Uc = evpHc* /c, hb^s the Bohr magneton, cr = ±1 
stands for spin up and down, respectively. 

It is important that Eq.(4) can be analytically solved: 



X exp 



exp[±iSex/va:{py)] 



±i 



HbO'Hx 



Vx{Pv) 



. Vx{Py) Jo 



cos I PzC* H u \ du 

Vf 



(6) 



where wave functions (6) are normalized on 5(ei — €2), 
Se = e—CF- The corresponding finite temperatures Green 
functions can be derived from Eq.(6) by means of the 



standard procedure [33]: 

9iijAx,Xl,Py,Pz,CT) = 



. sgn{ix>n) 
-I -, — r- exp 



UJn{x - Xi) 



MPy) 



X exp 
X exp 
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. Ux{Py) Jxi V VF 



(7) 



[Note that, in contrast to the previous works [26- 
30,11,20], Eqs.(6),(7) take into account dependence of 
electron velocity along conducting chains, Vx{py), on a 
momentum component Py. As shown below, it allows to 
describe accurately the Pauli paramagnetic destructive 
effects against superconductivity.] 

In this Letter, we consider a singlet d-wave like sce- 
nario of superconductivity in (TMTSF)2C104 conductor 
[14.17-19], which is consistent with all available experi- 
mental data. Therefore, we introduce the following su- 
perconducting order parameter. 



A{py,x) = V2cos{pyb')A{a 



(8) 



where the first term, y/2cos{pyb'), is responsible for the 
existence of zeros on QID FS, whereas the second term 
describes both the orbital effects against superconduc- 
tivity and possible LOFF like phase formation. Below, 
we derive a so-called gap equation for the superconduct- 
ing order parameter (8). using the Green functions (7). 
It is derived by means of the Gor'kov equations [33] for 
non-uniform superconductivity (see, for example, Refs. 
[34-36]. As a result of rather lengthly calculations, we 
obtain: 



A(.)=,'/4^/ 

J Vx(Py) J\x-xi 

J f StcVF 
[cJcVxiPy] 



2'!rTdxi 
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sin 
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2vf 
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MPy) 



2vf 

A{xi), 



(9) 



where g' stands for electron coupling constant, $7 is a 
cutoff energy, parameter /3 takes into account possible 
deviation of the so-called electron gi-factor, g = 2(3, from 
the value g = 2 [37]. We stress that Eq.(9) is different 
from the gap equations, used so far, and, unlike Refs. 
[26-30,11,20], it describes accurately not only the orbital 
effects but also the Pauli paramagnetic ones. Note that 
Eq.(9) is based on a quantum mechanical treatment of 
electron motion both in parallel and perpendicular to 
conducting layers directions. It is the most general gap 
equation, which can be written for QID conductor (2) 
in a magnetic field (3). It is possible to show that the 
major quantum parameter in Eq.(9) is 2tcVF/tOcVx{py) — 
2tc/u)c- It is also possible to prove that in low magnetic 
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fields, where 2tc/wc > 1 and (T^ - T)/Tc < 1, Eq.(9) is 
reduced to the well known GL equation [38]. 

Let us estimate a value of the dimensionless quantum 
parameter 1±{H) = 2tc/uJc in Eq.(9), which, using classi- 
cal language, represents a size of electron trajectory along 
z axis in terms of interlayer distance [26] : 



z{t, H) = lj_{H) c cos(wct), 
where t is time. It is easy to show that 

2\/2 tc 2 X 10^ tc tb 



TT aC*Hta H{T) tbta' 



(10) 



(11) 



where H{T) is a magnetic field, measured in Teslas. 
Here, according to Ref.[32], ta/h = 10 and, accord- 
ing to Ref.[38], tjh = {h*ly/2c*){H^JH^^)GL with 
{Hcil ^c2)gl is being a ratio of the GL slopes of the 
upper critical fields along c* and h' axes, correspond- 
ingly [39]. Note that the ratios ta/th = 10 [31] and 
i^cs/ ^c2)gl = 26 [12,13] are very well measured in 
(TMTSF)2C104 conductor. If we take H{T) = 6 T, we 
obtain 



1±{H = 6T) ~ 0.48, 



(12) 



which means that a size of electron classical trajectory 
along c* axis (10) is significantly less than interlayer 
distance, c* . In this case, which corresponds to the 
3D —5- 2D dimensional crossover of electron motion in 
a magnetic field [26,40], it is possible to make sure di- 
rectly from Eq.(9) that we can approximate the Bessel 
function as Jo{z) ~ 1 — 2^/4. 

Let us consider the above mentioned approximation for 
integral equation (9) at zero temperature, T = 0. It is 
possible to show that solution for a superconducting gap, 
A (a:), in this case can be written as 

A(a;) = exp(ifca;)[l-t-ai cos(2wca;/i'i;')-t-a2 sin(2wca;/w_F)], 

(13) 

where |ai|, |q:2| ^ 1- Eq.(9), determining the upper crit- 
ical field, in the same approximation and at T = can 
be expressed as 



1 
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, / UJcZ \ 

\2^J 



cos 



Vf 



-kz 



(14) 



where g is renormalized electron coupling constant, xi — 

X = zvx{py) /vf ■ [Note that we set ai = 0:2 = in 
Eq.(14), since we disregard all contributions of the order 
of 1\{H) < 1\{H) to the upper critical field.] 

Below, we simplify Eq.(14), taking into account that 
electron velocity component along conducting x axis is 



Vx{Py) = Vf[1 -)- acos(py6')], 



(15) 



where a = ^plthjta — 0.14 [20]. More specifically, 
Eq.(14) for a -C 1 can be written as follows: 
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' dz 
z 



2P^ibHz 



VF 



cos{kz)[Jo{akz) — J2{oLkz)\ 



X 



l-2/i(if)sin^(^ 



(16) 



It is important that Eq.(16) accurately takes into account 
the Pauli paramagnetic effects against superconductivity, 
unlike Refs. [26-30,11,20]. Note that, in the absence of the 
Pauli paramagnetic effects (i.e., at /3 = 0), Eq.(16) de- 
scribes the Recintrant superconducting phase [26] with 
dTc/dH > 0. Therefore, we call superconducting phase, 
described by Eqs.(16),(18), the hidden Reentrant super- 
conductivity. 

Let us further simplify Eq.(16) by taking into account 
that 



2'KTrdz 



(17) 



where T,. is the superconducting transition temperature 
at H = 0. As a result, we obtain 



x< cos{kz)[Jo{akz) — J2{akz)] 



l-2fj_{H)shr 



\2vfJ 



(18) 



where fisH* = ttTc/2j, 7 is the Euler constant. Numer- 
ical analysis of Eq.(18) shows that the upper critical field 
along b' axis, H^2, for 1±{H) = 0.48 and /3 = 0.84 has a 
maximum at = 0.88(2/3/isiJ/wir) and is equal to 



H^^2 =^ 5.9 T. 



(19) 



[We pay attention that the obtained value of the upper 
critical field (19) well corresponds to the value of a mag- 
netic field (12).] For the same values of the parameters 
1±{H) and /?, numerical analysis of Eq.(9) gives the fol- 
lowing values for factors ai and 02 in Eq.(13): 



ai 



-0.139, Q!2 = 0.021 i. 



(20) 



Below, we summarize the main results of the Letter. 
We have derived gap equations (9), (14), (16), (18), which, 
unlike gap equations in the previous publications, take 
accurately into account not only the orbital effects, but 
also the Pauli paramagnetic effects against superconduc- 
tivity. We have analyzed the experimental data [12,13,32] 
and shown that, in contrast to the common believe, the 
quantum effects of electron motion in a magnetic field 
[26,41] are strong in relatively weak magnetic fields of 
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the order of 5 - 6 T in (TMTSF)2C104 conductor. By 
analyzing the above mentioned gap equations, we have 
explained how superconductivity in (TMTSF)2C104 can 
exceed both the quasi-classical upper critical field [24] 
and Clogston paramagnetic limit [25] and how it can 
reach its experimental value, if ~ 6 T [12,13]. We 
have shown that, due to the reentrant quantum effects 
[26,41], superconductivity survives in form of the hid- 
den Reentrant superconducting phase, corresponding to 
thrc;c LOFF-likc phases. Although we have not calcu- 
lated in the Letter phase diagram of the (TMTSF)2C104 
superconductor in all range of temperatures and mag- 
netic fields, we anticipate the existence of phase transi- 
tion between the BCS and LOFF phases at iJ ~ 2.5 T, 
which can be experimentally studied. 

Ill conclusion, we note that the considered above hid- 
den Reentrant superconductivity is a rather general phe- 
nomenon. It is expected to exist in other (TMTSF)2X 
conductors and may exist in quasi-two-dimensional su- 
perconductors in a parallel magnetic field. Neverthe- 
less, this phase in (TMTSF)2PF6 material, which is sta- 
ble in a mixed superconducting-spin-dcnsity-wave state 
[21,22,42] in a magnetic field up to = 9 T, possesses 
some peculiarities. Our preliminary analysis shows that, 
to describe the hidden Reentrant superconducting phase 
in (TMTSF)2PF6, it is necessary to take into account 
some additional effects such as the singlet-triplet mixing 
phenomenon [43] or possible singlet-triplet phase transi- 
tion (see, for example, [30,20]). 
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